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The  optimal  linear  quadratic  guidance  problem  with  an  accelerating 
target  is  investigated.  The  missile  state  is  partitioned  into  a kine- 
matic state  and  an  airframe  state.  Both  the  penalty-weighted  and  con- 
strained terminal  state  cases  are  treated.  The  resulting  optimal  guid- 
ance law  requires  estimates  of  target  accelerations  which  are  derived 
via  linear  observers.  Results  of  a point  mass  missile,  a one-time 
constant  missile,  and  a two-time  constant  missile  are  given. 
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SECTION  I 


INTRODUCTION 


This  report  examines  minimum  control  cost,  minimum  terminal  miss 
guidance  for  the  intercept  of  a maneuvering  target  by  a missile  with 
inherent  airframe  and  autopilot  dynamic  properties.  Particular  atten- 
tion is  given  to  the  idealized  problem  of  zero  terminal  miss,  wherein 
the  control  gains  are  given  in  terms  of  the  state  transition  of  the  un- 
coupled airframe  dynamics.  This  approach  separates  the  kinematic  por- 
tion of  the  intercept  dynamics,  which  is  common  for  all  intercept  prob- 
lems, from  the  kinetic  portion  of  the  missile  dynamics.  The  particular 
problem  structure  makes  the  results  applicable  to  a variety  of  pursuit- 
evasion  problems  in  which  the  missile  airframe  can  be  represented  by  a 
linear  model. 

The  resulting  control  law  developed  in  this  paper  has  a term  re- 
lated to  the  intercept  kinematics,  which  is  recognizable  as  the  general- 
ized proportional  navigation  term.  A second  term  in  the  control  law  is 
a linear  function  of  the  missile  airframe  state  and  represents  the  guid- 
ance compensation  due  to  finite  airframe  dynamics.  The  final  term  in 
the  guidance  law  is  related  to  target  motion,  providing  an  effective 
control  in  cases  where  target  motion  can  be  measured  or  predicted 
accurately. 

The  next  section  outlines  the  problem.  Subsequent  sections  develop 
a minimum  miss,  and  a zero  miss  control  for  the  general  problem.  A 
later  section  treats  the  estimation  of  the  unknown  target  acceleration. 
Finally,  the  analysis  is  carried  to  conclusion  for  the  case  where  the 

kinetics  of  the  airframe  are  independent  of  the  kinematic  state. 
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SECTION  II 


FORMULATION  OF  THE  GUIDANCE  PROBLEM 


This  formulation  structures  the  guidance  problem  to  separate  the 
kinematics  of  the  intercept  from  the  dynamics  of  the  missile.  The  motion 
of  the  target  is  accepted  as  an  uncontrollable  input  to  the  problem;  how- 
ever, the  kinetic  state  equation  can  be  augmented  with  a target  model. 
Consider  first  the  intercept  kinematics. 


1.  KINEMATICS 

Let  the  vector  position  and  velocity  of  the  target  (T)  and  inter- 
ceptor missile  (I)  be  represented  in  an  inertial  frame  by  y and  v,  as 
illustrated  in  Figure  1. 


y 

V 


I 

I 


V * V 

yT  T 


VT  ~ aT 


(1) 


Defining  the  relative  position  and  velocity  of  the  target  with  respect 
to  the  missile  yields 


y = v 

v = aT  - aT 


(2) 


Letting  the  vector  x represent  the  kinematic  state  of  the  intercept, 
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then 

x - Ax  + B (ax  - aT)  (3) 

where 


In  (3)  the  identity  and  null  submatrices  reflect  the  dimension  of  the 
problem. 

2.  KINETICS 

The  airframe/autopilot  response  state  is  designated  as  z and  satis- 
fies the  linear  equation 


z = Dx  + Ez  + Fu  (4) 

subject  to  the  airframe  control  u (thrust,  control  surface  deflection, 
etc.).  The  dimensions  and  components  of  the  coefficient  matrices  in 
(4)  are  missile  dependent  and  provide  the  generality  in  the  problem. 

The  intercept  kinematics  are  coupled  to  the  airframe  dynamics  by 


ax  = Gz  + Hu  (5) 

; 

$ Thus  any  linearized  airframe  describable  by  (4)  and  (5)  is  subject  to 

,,  the  analysis. 

% 

3.  OPTIMAL  GUIDANCE 

The  conventional  guidance  performance  index  for  homing  missiles  is 
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cost  and  time  respectively.  In  those  cases  where  the  terminal  miss  is 
the  significant  cost,  it  is  logical  to  attempt  to  constrain  the  final 


position  y(t^)  to  zero  and  to  develop  the  corresponding  control  law 
under  this  condition.  Equation  (6)  may  be  replaced  by 


Qsu  Ru  + y)dt 


Tx(tf ) = 0 


-b  o] 


The  equations  (3-7)  are  collected  below. 


Minimum  Miss 


x = Ax  + BGz  + BHu  - Ba„ 


z = Dx  + Ez  + Fu 


J = Jsx(tf)T  Sfx(tf)  + I (*suT  Ru  + y)dt 


Zero  Miss 


x = Ax  + BGz  + BHu  - Be 


z = Dx  + Ez  + Fu 


Tx(tf ) = 0 


(*SU  Ru  + y)dt 


SECTION  III 


MINIMUM  TERMINAL  MISS 


If  the  performance  index  in  (8)  is  augmented  in  the  usual  manner, 
the  Hamiltonian  is 


H = *suT  Ru  + y + XT  [Ax  + BGz  + BHu  - Ba^]  + p [Dx  + Ez  + Fu]  (10) 


The  resulting  boundary  value  problem  is 


x = Ax  + BGz  + BHu  - Ba„ 


x(to) 


z = Dx  + Ez  + Fu 


* T T 

X = -A  X - D M 


• m m rp 

M = - G B x - E b 
u = -R-1  (HTBT  X + FTm) 
H(tf)  = 0 


z(t0) 


X(tf) 


w(tf) 


sf  x(tf) 


The  computation  of  the  control  gains  is  achieved  via  the  inverse 


formulation 


x = Qx  X + Q2  p + Q3 


z = X + Q5  y + Q6 


leading  to  the  equations 


• T T T T —1  T T 

Q1  = AQX  + QjA  + Q2G  B + BGQ2  - BHR  H B 


Ql(tf)  = sf- 


Q2  = AQ2  + Q2ET  + Q1DT  + BGQ5  - BHR-1  FT 
Q3  = AQ3  + BGQ6  - BaT 


Q2(tf)  = 0 


Q*-Q2 


Q5  - EQ5  + Q5ET  + DQ2  + Q2TDT  - FR-1  FT 


Q5(tf)  - sf 


(12) 

(Contd. ) 
*-l 


q6  - EQ6  + DQ3 


Q6(t£)  - 0 


The  nonsingular  diagonal  matrix  Is  used  for  the  computation  of 
the  inverse  problem  and  its  elements  are  set  to  zero  in  the  solution  for 
the  gains. 

The  resulting  control  vector  u is 


u = + K2z  + K3 


where 


and 


— ITT  T T 

K1  = -R  (H  B P1  + F P2  ) 


-ITT  T 

K2  = -R  (H  B P2  + F P5) 


-ITT  T 

k3  = -R  (hVp3  + fap6) 


P1  ■ “>!  - ^s'1  ,>2Tl'1 


(13) 


p2  • -w1  - OjQ,-1  a2Vx  82fi5‘x 


P5  - lq5  - q2T  h1  q21'1 


(U) 


P3  " - P1Q3  ' Vfi 


P6  " " P2  Q3  “ P5Q6 


The  particular  case  of  interest  is  that  in  which  the  kinetics  of 
the  airframe  are  indeper  ’"ct  of  the  intercept  position  and  velocity,  i.e. 


L 


f 


1 


L 


D = 0. 

In  Lnis  case,  the  equations  (12)  are  easily  integrable  yielding 


Qs 


/f 

I $E  (t-x)  FR_1  FT4i. 

J: 


eT  (t-T)  dT  + S*"1 


A (t-T)  B [HR-1  FT  - GQ5  (t)]  *ET(t-T)  dx 


Q1  = M A (t_T)  [BHR_1  hT®T  " BGQ2T  " T (f-"T)  dT  + s 


r 


(15) 

-1 


Q6  = 0 


1} 


Q3  - - I *A  ^t-T^  BaT  ^T^dT 


where  and  are  the  transition  matrices  for  A.  and  E. 


*A  = 


I tl 
0 I 


and 

= E»E  » ^E(0)  = I* 

The  elements  of  S*  are  set  to  zero  in  the  resulting  expressions 
for  the  elements  of  Q. 

The  minimum  terminal  miss  control  law  is  determined  by  (13)  - (15) 
The  case  where  the  terminal  miss  weighting  is  large  is  solved  by  con- 
straining the  terminal  position,  as  is  done  in  the  next  section. 


* JB*  , 


SECTION  IV 


ZERO  TERMINAL* MISS 

The  zero  terminal  miss,  minimum  control  coBt,  minimum  time  problem 
formulated  in  an  earlier  section  and  outlined  in  (9)  is  solved  here. 
This  problem  is  also  treated  m reference  1.  The  augmented  performance 
index  for  the  problem 

te 


J = vT  Tx(t^)  + ^ lkuT  Ru  + y + XT  (Ax  + BGz  + BHu  - BaT) 


+ p (Dx  + Ez  + Fu)]dt 

yields  the  boundary  value  problem  described  by 


Ax  + BGz  + BHu  - BaT 

x(to) 

= X 

o 

Dx  + Ez  + Fu 

z(to) 

= z 

o 

-atx  - dt 

*(tf) 

T 

= Tv 

-gtbtx  - etu 

M(tf) 

= 0 

-r_1(htbtx  + FTp) 

H(tf  ) 

= 0 

(16) 


Selecting 


(16)  becomes 


X = S, x + S„z  + S0v  + S. 
1 2 3 A 


h = S^x  + Sfiz  + S^v  + Sg 


r. 


1 ,-13— a.. 


«.  * 4 .44< 


Inspection  of  (17)  and  (18)  shows  that  all  unknown  matrices  except 
and  are  null.  Leaving 


• T T 

S3  + A S3  + D S?  = 0 


• T T T 

S?  + E S?  + G B S = 0 


- 1 T T T 

u = -R  (H  B S3  + F S?)v 


The  invariance  of  the  terminal  manifold 


S3(tf)  - T 


S?(tf)  = 0 


Y(t)  = S9X  + S1QZ  + S13v  + S12 


implies 


S9  + S9A  + S10D  = 0 


S10  + S10E  + S9BG  = ° 


S9(tf)  = T 


S10(tf>  = ° 


Sn  “ (SgBH  + S1qF)  R 1 (HTBTS3  + FTS?)  = 0 Sn(tf)  = 0 


B12  ~ B9BaT  = G 


S12(tf>  ' 0 


where 


V = "Sll  [S9X  + S10Z  + S12] 


Therefore  the  control  law  for  zero  miss  is 


u = R_1  (HTBTS3  + FTS?)  Sn_1  (S9x  + S1Qz  + S12), 


where  the  navigation,  airframe,  and  target  motion  components  are  evident. 


r 6 ■'  . . ' : " 
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For  the  uncoupled  airframe  case  where  D - 0,  the  equations  (19)  and 
(20)  are  integrated  for  computation  of  the  optimal  gains.  Integration 
yields  f-  T 1 


(tf  -t)I 


S (t)  BG4>*  (t-T)  dT 
E 


-tf  (21) 

- J lS1Q(r)  F + S9(t)  BH]  R_1  [S1q(t)  F + S9(t)  BH]T  dr 


* * 

where  $ = -$  E 

£ £ 


The  target  maneuvering  term  of  the  control  for  the  uncoupled  case 


can  be  written 


where 


ms.  ■*  m 


Also 


(23) 


K3  = K30  (tf  ' ° tVT(t)  " VT(t)1 

where  VT(t)  represents  present  target  velocity  and  V(t)  Is  the  average 
relative  velocity  on  the  terminal  interval  (t,  tf) . The  need  for 
effective  estimation  of  target  motion  is  recognized. 
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SECTION  V 

TARGET  MANEUVERING 

For  minimum  miss  control  the  target  maneuver  term  of  the  control  is 
uT  - -KjQj  - K2Q6  (24> 


Kx  = -R_1(HTBTP1+FTP2T) 


— 1 T T T 
K2  = -R  1(HiBiP2+F  P5) 


Q3  = %33(t-T)  BaT(T)dT 


Q6  = I $63(t-T)  BaT(T)dr 


33 

36 

63 

% 

66 

is  the  state  transition  matrix  corresponding  to 


A BG 


D E 


For  zero  miss  guidance  the  target  maneuver  component  of  the 
control  reduces  to 

f£ 

u^,  = -R_1(HTBTS3+FTS7)  I S9(T)BaT(r)dT. 

14  ^ 


' 

_ffl r1"  -v  1 -■  --1— - 


j 


In  both  the  minimum  miss  and  zero  miss  problems,  the  segment  of  the 
control  due  to  the  maneuvering  target  requires  a weighted  integral  of 
aT>  This  control  can  be  implemented  for  the  deterministic  case  if  a 
reliable  estimate  of  aT  can  be  obtained  from  an  on-line  reconstructor. 
Such  an  estimate  may  be  taken  from  a linear  observer,  if  the  unknown 
target  acceleration  is  modeled  by  a linear  equation 


aT  = Mi3t  + M2 


This  assumption  is  partially  justified  by  consideration  of  the  dominant 
time  constants  of  the  target  compared  to  the  normal  intercept  time.  In 
most  cases  a constant  or  ramp  acceleration  adequately  describes  the 
target  motion  in  this  interval.  If  the  states  x and  z are  completely 
measurab le. 


a simple  reduced-order  linear  observer  for  a^  may  be  used.  To  facili- 
tate the  development,  this  approach  is  taken,  although  a higher  order 
observer  can  be  used  if  components  of  x and  z are  not  directly  measur- 
able but  are  observable. 

The  equations  are 


x = Ax  + BGZ  + BHu  - Ba„ 


z = Dx  + Ez  + Fu 


aT  = MxaT  + M2 


If  an  estimate  a^  satisfies 


a^  = a„  + Lx 
T T 


15 


then 

i 

aT  - a^.  = (LB+M^)  (29) 

and  the  matrix  L may  be  chosen  to  give  the  desired  convergence  for  the 
observer.  In  the  actual  implementation,  line-of-sight  angles  and  rates 


will  replace  position  and  velocity  in  the  observer  model,  and  some  com- 


ponents of  z can  be  added  to  the  observed  state,  as  mentioned. 


I 


SECTION  VI 


TARGET  INTERCEPT  SITUATIONS 
USING 

OPTIMAL  GUIDANCE  WITH  TARGET  ESTIMATION 


Situation  I Kinematic  Example 

A preliminary  example  is  one  of  a longitudinal  planar  intercept  of 
a target  with  an  unknown,  constant  acceleration  by  a point  mass  missile. 
To  facilitate  the  solution,  a zero  miss  constraint  with  quadratic 
missile  acceleration  cost  is  assumed.  Since  the  lateral  motion  is  com- 
pletely uncoupled,  the  same  type  of  control  is  applicable  in  that 


plane . 


The  equations  are 


y = v 


J = h 


o 1 


0 0 


, y(tf)  = 0 


D = 0 


H = 1 


R = 1 


Integration  yields  from  (21) 


S9  = C°  Sri 


S11  = -1/3  V 


S 12  = I (tf-r)  aT(r)dT 


yCi-»4-  r.-.  TH 


where  t ■ t,  - t 

g f 


The  control  is 


aI  = (3/tK2>y  + (3/Ov  + (3/t  2)  I (tf-t)  ax(T)dT 


If  the  coordinate  axes  are  chosen  such  that  the  down-range  coordinate 
is  along  a non-rotating  line  of  sight,  the  line-of-sight  angle  A is 


A %y(t)/rg(t) 


where  r^(t)  Is  the  range-to-go,  which  depends  on  the  target  and 
missile  accelerations,  as  well  as  present  relative  velocity  and 
position. 


y = Ar 


v = Ar  + Ar 
g 8 

Hence  t^ 

aI  = (3/tg2)(rg+^gtg)  + (3rg7tg)X  + (3/tg2)  f (tf'T)  aT(T)dT 


If  the  range-to-go  is  decreasing  steadily  along  the  line  of  sight, 


r = r t 
g o g 


the  commanded  acceleration  becomes 


aT  = 3r  A + (3/t  2)  | (t,-r)  aT(t)dr 

1 o g f T 

In  any  case,  if  the  control  for  intercept  of  a maneuvering  target  is 
required,  the  target  acceleration  normal  to  the  line  of  sight  a^.  must 
he  estimated.  If  a^  is  an  unknown  constant  as  specified,  it  is 


modeled 


Hence 


aT  = 0 


= 0 


“ 0 


In  this  case  the  target  observer  model  is 


aT  = aT  + *2v 


~^2ax  ^2al 


where  l„  may  be  taken  large  enough  to  satisfy  convergence  requirements 


for  the  observer. 


If  a constant  estimate  of  a^  is  achieved,  then 


= (3/t  2)y  + (3/t  )v  + 3/2aT 


Due  to  the  terminal  constraint,  the  miss  distance  will  be  reduced  to 

zero  no  matter  what  error  exists  in  the  target  estimate  (at,  perhaps, 

a prohibitive  terminal  control  gain).  The  control  cost  will,  however, 

reflect  the  error  in  a . 

T 

For  comparison,  assume  that 


aT  = Y3T 


a:  - (3/t  2)y  + (3/t  )v  + (3/2)yaT 


Thus  y = 0 corresponds  to  ignoring  target  maneuvers,  and  y = 1 cor- 
responds to  optimal  control.  The  reduction  in  control  cost  from  that 


* i"i  i iBh I iTT' ii ' 1 m in  ■ . • '..***»>  . f.. . 


obtained  by  not  using  a target  estimate  is 


AJ  = 9(2-y)YaT2tf  /8 


While  any  target  estimate  (for  constant  accelerating  target)  such 
that  0<y<2  represents  an  improvement  in  cost,  the  minimum  cost  is 
obviously  achieved  when  the  estimate  is  perfect. 

V 

Situation  II  One-Time-Constant  Missile 

A more  realistic  situation  occurs  when  an  airframe  state  exists  in 
the  missile  model.  In  this  case  the  equations  of  intercept  are 

y = v 
v = aT  - 
z = a(u  - z) 

az  = gz 


Referring  co  (3,  4,  5) 


r—  — 

0 1 

0 

A = 

B = 

_o  o_ 

-1 

D = (0  0)  E = -a 

F = a G = g 

H = 0 

Setting 

R = X 

and 

s = 1 
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The  substitution  yields 

u = K^y  + K^v  + K2z  + K3 


where 


-2  -1  -2  atg 

K„  = K (a  - a t - a e °) 

2 11  g 


K3  - Kn  I (tf-T)  aT(i)dx 


If  a reduced  observer  is  used  to  reconstruct  the  target  acceleration 
for  equations  (27,  28)  yield 


aT  = aT  + *ly  + i2V 


aT  = (-)l2+ml^aT  “ 5'1V  + ^2gZ  + m2 


where  the  error  e = a^  - a^,  decays  according  to 


e = ~®'2e 


If  the  target  acceleration  is  almost  constant  over  the  duration  of  the 
intercept,  and  are  taken  as  zero.  If  Z,  is  taken  to  be  zero  for 


convenience , 


aT  = aT  + V 


= "*2^  + J'2gZ 


The  observer  is  shown  in  Figure  2. 
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Figure  3.  Two-Time-Cons  tan c Missile 
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If  the  initial  line  of  sight  is  taken  as  the  x-axis,  for  small  9 
the  equations  of  intercept  are  as  follows: 


yT  • y!  = y 


VTx  VIx 


v = 0 
o 


V — V = v 

Ty  VIy 


v = (T/M) 6 
j-y 


..  v.  , c 2 

9 = (7^)9  + ( ■)  u 
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m 


yielding 


y = v 


v = g0  + a 


0 = 0) 


o)  = aw  + Bu 


An  optimal  guidance  problem  is  formulated  with 


J = y(tf)2  + A I u2dt 


In  the  format  of  Section  II, 


G = (-g  0) 


H = 0 


Solving  for  the  Q and  subsequent  gain  matrices  in  (13-15) 


W ■ "u'V  [1 


K2(tg>  - K11<tg)  [-  “'«] 


- K„(t  > (tf-t)  ^(OdT 


Figures  13-20 


give  typical  trajectories  for  this  intercept 
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Figure  5.  Miss  Distance  for  One-Time  Constant  Missile. 

Target  Acceleration  = 2g 
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Figure  6.  Miss  Distance  for  One-Time  Constant  Missile. 
Target  Acceleration  = 4g 
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Figure  9.  Miss  Distance  for  One-Time  Constant  Missile 
Target  Acceleration  = -4g 
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Figure  13.  Miss  Distance  for  Two-Time  Constant  Missile 
Without  Target  Estimation 
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Figure  15.  Miss  Distance  for  Two-Time  Constant  Missile. 
Target  Acceleration  = 4g 
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f Figure  17.  Miss  Distance  for  Two-Time  Constant  Missile 

Target  Acceleration  = -2g 
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Figure  19.  Miss  Distance  for  Two-Time  Constant  Missile. 
Target  Acceleration  = -6g 
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Figure  20.  Terminal  Miss  Distance  for  Two-Time  Constant  Missile 
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